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The  Effects  of  Insulator  Wall  Material  on  Hall  Thruster  Discharges: 

A  Numerical  Study  ' 


John  Michael  Fife 
Air  Force  Research  Laboratory 
Edwards  AFB.  California 

Summer  Locke 
University  of  Washington 
Seattle.  Washington 


An  investigation  was  untiertaken  to  determine  how  the  choice  of  insulator  wall  material  inside  a  Hall  thruster 

discharge  channel  might  affect  thruster  operation.  In  order  to  study  this,  an  evolved  hybrid  particle-in-cel!  (PIC) 

numerical  Hall  thruster  model.  HPHall,  was  used.  HPHall  solves  a  set  of  quasi-one-dimensional  fluid  equations  for 

electrons,  and  tracks  heavy  particles  using  a  PIC  method.  The  hvo  systems  are  linked  by  charge  neutrality.  Five  cases 

were  executed  with  various  secondary  electron  emission  coefficients  at  the  insulator  walk  Results  show  a  steady  .'VVirAs 

increase  in  thruster  efficiency  from  0.434  to  0.483  as  secondary  electron  emission  coefficient  is  decreased  by  28.  This 

suggests  that  secondary  electron  emission  be  one  of  the  parameters  considered  during  the  Hall  thruster  design  phase, 

and  that  channel  insulator  materials  with  low  secondary  electron  emission  should  be  favored. 


Introduction 


Electrons 


The  choice  of  ceramic  insulator  material  for  Hall 
thruster  discharge  chambers  is  largely  dependent  on 
structural  and  thermal  requirements,  as  well  as 
resistance  to  ion  sputtering.  Analysis  has  recently 
indicated  that  the  secondary  electron  emission 
coefficient  of  the  discharge  channel  insulator  may 
affect  the  plasma  discharge  and.  therefore,  also  be  an 
important  Hall  thruster  design  parameter.  This  paper 
investigates  changes  in  the  Hail  thruster  discharge  for 
various  secondary  electron  emission  coefficients  using 
a  numerical  Hall  thruster  simulation,  HPHall. 

Several  numerical  simulations  of  Hall  thrusters  have 
been  developed  by  other  reseachers.  1-2  HPHall  is  a 
two-dimensional  transient  hybrid  paiticle-in-cell 
(hybrid  PIC)  simulation.  Although  descriptions  of 
HPHall  have  been  presented  before, w  X  brief 
summary  of  the  model  and  method  will  be  repeated 
here. 

Governing  Equations 

Although  this  transient  2-D  3-V  simulation  operates  in 
cylindrical  coordinates  ( z.r.v..vr,ve ),  some  vector 
quantities  in  the  analysis  below  are  written  with 
res(iect  to  the  magnetic  Field  lines.  As  Fig.  1  shows. 
h  and  /  are  used  to  represent  the  distance  vectors 
normal  and  tangent  to  the  magnetic  field  lines, 
respectively. 


The  simplified  electron  equations  consist  of  a 
generalized  Ohm’s  law,  a  current  conservation 
equation,  and  an  electron  temperature  equation. 
Assuming  a  Maxwellian  electron  distribution, 
quasineutrality,  and  a  particular  ion  field,  these  three 
equations  are  sufficient  to  yield  electron  current 
density,  space  potential,  and  electron  temperature  as  a 
function  of  time. 

The  diffusion  coefficient  of  electrons  along  magnetic 
field  lines  is  assumed  to  be  much  greater  than  the 
diffusion  coefficient  across  them.  Ignoring  the 
magnetic  mirror  effect,  and  assuming  constant 
electron  temperature  along  magnetic  field  lines,  the 
momentum  equation  gives, 

(I) 

Note  that  Eq.  (1)  holds  along  magnetic  field  lines.  A 
is  a  magnetic  stream  function  which  is  constant  for 
any  given  field  line. 


Electron  diffusion  across  the  magnetic  field  is 
assumed  to  obey  a  Generalized  Ohm’s  Law.  In  the  lab 
frame,  the  cross-field  electron  velocity  in  terms  of  an 
effective  electron  mobility  across  the  magnetic  field 
lines.  ft'±,\ s 


(2) 


Fig.  1.  Diagram  showing  a  simplified  Hall  thruster 
discharge  channel  cross  section,  the  coordinate  system 
used,  and  the  current  convention. 

Classical  cross-field  mobility  in  the  weakly  ionized 
limit7  (~l/fi2)  has  been  shown*  not  to  adequately 
describe  the  high  electron  transport  across  lines  of 
force  in  the  presence  of  a  strong  magnetic  field.  Some 
previous  work  has  suggested  that  the  discrepancy 
between  measured  and  predicted  mobility  may  be  due 
to  anomalous  “Bohm”  diffusion.9  which  goes  as  MB . 
Neither  classical  nor  Bohm  models  appropriately 
describe  the  mobility  throughout  the  discharge 
channel^Jt  is  believed  thafTiin  some  regions,  other 
"mechanisms^  such  as  walK conductivity  or  wave 
transportf'may  dominate.  However,  in  the  interest  of 
exploring  the  physics  of  the  acceleration  process,  we 
use  classical  and  Bohm  mobility  here. 


(3) 


and  leave  the  coefficient.  K,,  as  an  adjustable 
parameter  between  0  and  1.  Based  on  comparison 
with  experimental  data  the  “best”  K,  was  found  to 
be. 15. 


Current  Conservation 

Since  quasineutrality  is  imposed,  no  space  charge  can 
accumulate,  and  current  must  be  conserved  for  the 
whole  device.  A  conservation  equation  for  current 
crossing  any  magnetic  field  line  can  be  written 
l„  ~  I,  +  + /,  .  where  /, ,  and  /.  are  the 

discharge,  electron,  ion.  and  near-wail  currents, 
respectively.  In  terms  of  integrals  along  magnetic 
field  lines. 


By  combining  Eqs.  (1)  through  (4)  a  differential 
equation  may  be  obtained  for  the  cross-field  variation 
of**. 

Electron  Energy  Equation 

An  electron  energy  equation  is  derived  under  the 
assumptions  that  electrons  have  a  Maxwellian  velocity 
distribution,  and  that  the  pressure  dyad  reduces  to  a 
scalar  pressure  term;  nrkTt.  Source  terms  include 

losses  due  to  ionization,  radiation,  and  charge-field 
interactions.  Ionization  and  radiation  losses  are 
modeled  analytically  with  a  net  ion  production  cost 
according  to  Dugan  and  Sovie.10  Charge-field 
interactions  are  modeled  as  j)lc,.  The  energy 

equation  is  applied  across  magnetic  field  lines  only. 
Along  them,  electrons  are  assumed  to  be  isothermal. 

Ionization  Rate 

The  bulk  electron-neutral  ionization  rate  is  determined 
by  integrating  the  Drawin7  cross-section  over  a 
Maxwellian  electron  distribution.  This  bulk  ionization 
rate,  plus  the  equations  of  motion  for  the  ions  and 
neutrals,  completes  the  model  for  those  species. 

Modeling  Wall  Effects 

Interaction  of  the  Hail  thruster  plasma  with  the 
insulator  wall  depends  upon  secondary  electron 
emission  coefficient  of  the  wall  material.  Previous 
models3  assumed  zero  or  constant  secondary  electron 
emission  coefficient,  resulting  in  very  low  electron 
energy  flux  to  the  wall.  This  caused  high  predictions 
of  T, .  The  following  theory  is  a  more  detailed 
approach,  parts  of  which  have  been  published 
before.5-6-11  Also.  Katz  et  al.12  present  an  analogous 
theory  for  computation  of  charging  rate  on 
geosynchronous  spacecraft. 

Experimental  data  show  that  the  ratio  of  secondary  to 
primary  electrons,  S ,  is  a  function  of  incident  electron 
energy,  £  (eV).  At  low  (ciOOeV)  energies,  an 
exponential  fit  can  be  used  to  closely  approximate 

S(E): 

S-AE'  (5) 

first  we  look  at  a  negative  wall  potential  with  respect 
to  the  plasma  (*„<0),  and  an  ion-attracting  sheath. 
From  sheath  theory,  the  primary  electron  Dux  is. 


I>-2£,  (6) 

4 

where  c,  is  the  electron  mean  thermal  speed.  The 
integral. 
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Fig.  1.  Diagram  showing  a  simplified  Hall  thruster 
discharge  channel  cross  section,  the  coordinate  system 
used,  and  the  current  convention. 

Classical  cross-field  mobility  in  the  weakly  ionized 
limit1  (-1  IB2)  has  been  shown8  not  to  adequately 
describe  the  high  electron  transport  across  lines  of 
force  in  the  presence  of  a  strong  magnetic  field.  Some 
previous  work  has  suggested  that  the  discrepancy 
between  measured  and  predicted  mobility  may  be  due 
to  anomalous  “Bohm”  diffusion.9  which  goes  as  MB . 
Neither  classical  nor  Bohm  models  appropriately 
describe  the  mobility  throughout  the  discharge 
channel.  It  is  believed  that,  in  some  regions,  other 
mechanisms,  such  as  wall  conductivity  or  wave 
transport,  may  dominate.  However,  in  the  interest  of 
exploring  the  physics  of  the  acceleration  process,  we 
use  classical  and  Bohm  mobility  here. 

,3> 

and  leave  the  coefficient.  K ,,  as  an  adjustable 
parameter  between  0  and  I.  Based  on  comparison 
with  experimental  data,  the  “best”  K .  was  found  to 

Current  Conservation 

Since  quasineutrality  is  imposed,  no  space  charge  can 
accumulate,  and  current  must  be  conserved  for  the 
whole  device.  A  conservation  equation  for  current 
crossing  any  magnetic  field  line  can  be  written 
/,  =  I,  +  +  /„ ,  where  /, .  lr.  I, ,  and  /„  are  the 
discharge,  electron,  ion.  and  near-wall  aments, 
respectively.  In  terms  of  integrals  along  magnetic 
field  fines. 

i  i 

/„  =  -2zej  ntueJ  rds + 2 Kej  neu;JI  rds  +  Ir .  (4) 


By  combining  Eqs.  (1)  through  (4)  a  differential 
equation  may  be  obtained  for  the  cross-field  variation 
of  0*. 

Electron  Energy  Equation 

An  electron  energy  equation  is  derived  under  the 
assumptions  that  electrons  have.a  Maxwellian  velocity 
distribution,  and  that  the  pressure  dyad  reduces  to  a 
scalar  pressure  term.  ntkTt.  Source  terms  include 

losses  due  to  ionization,  radiation,  and  charge-field 
interactions.  Ionization  and  radiation  losses  are 
modeled  analytically  with  a  net  ion  production  cost 
according  to  Dugan  and  So  vie.1"  Chatge-field 
interactions  are  modeled  as  j]  I oc .  The  energy 
equation  is  applied  across  magnetic  field  lines  only. 
Along  them,  electrons  are  assumed  to  be  isothermal. 

Ionization  Rate 

The  bulk  electron-neutral  ionization  rate  is  determined 
by  integrating  the  Drawjn7  cross-section  over  a 
Maxwellian  electron  distribution.  This  bulk  ionization 
rate,  plus  the  equations  of  motion  for  the  ions  and 
neutrals,  completes  the  model'for  those  species. 

Modeling  Wall  Effects 

Interaction  of  the  Hall  thruster  plasma  with  the 
insulator  wail  depends  upon  secondary  electron 
emission  coefficient  of  the  wail  material.  Previous 
models3  assumed  zero  or  constant  secondary  electron 
emission  coefficient,  resulting  in  very  low  electron 
energy  flux  to  the  wall.  This  caused  high  predictions 
of  Tr .  The  following  theory  is  a  more  detailed 
approach,  parts  of  which  have  been  published 
before. 3-611  Also.  Katz  et  al.lj  present  an  analogous 
theory  for  computation  of  charging  rate  on 
geosynchronous  spacecraft. 

Experimental  data  show  that  the  ratio  of  secondary  to 
primary  electrons,  5 ,  is  a  function  of  incident  electron 
energy,  E  (eV).  At  low  (<l00eV)  energies,  an. 
exponential  fit  can  be  used  to  closely  approximate 

5(E): 

5  =  AE*  (5) 

First  we  look  at  a  negative  wall  potential  with  respect 
to  the  plasma  ( ^  <  0 ),  and  an  ion-attracting  sheath. 
From  sheath  theory,  the  primary  electron  flux  is. 


where  ~c,  is  the  electron  mean  thermal  speed.  The 
integral, 
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0  ■* 


(7) 

yields  the  secondary  electron  flux.  ,  and 

the  effective  secondary  emission  coefficient.  8# .  in 
the  form: 

8„]=T[2  +  B)AlkTtle)\  (8) 

where  F[.r]  is  Euler's  Gamma  function. 

imposing  wall  neutrality  by  balancing  the  fluxes  of 
ions  and  primary  and  secondary  electrons,  it  can  be 
shown  that 

<t=^in[(!-^)e,/3-S-],  (9) 

e  4i  ■„ 


0  10  T,-*,20  30  40 


Electron  Temperature  (eV) 

Fig.  2.  Normalized  heat  flux  to  the  wall  versus 
electron  temperature  for  xenon  on  boron  nitride. 


where  the  Bohm  velocity  is. 


From  (9).  the  sheath  potential  reverses  at  a  breakpoint 
temperature.  Telie, .  where. 

(11) 

The  wall  potential  is  negative  for  Tt<Tegf,  and 
positive  for  Tr  >  Tr  BP .  For  xenon  and  boron  nitride  (a 
candidate  insulator  for  Hall  thrusters,  for  which 
^  =0.141  and  £  =  0.576 ).B  T,„-  16.55eV.  and 
^.=0.997. 

For  the  case  of  the  ion-repelling  sheath,  the  neutrality 
condition  requires  that  8#  *  I .  Since  secondaries 
must  now  overcome  the  sheath  barrier  to  escape  into 
the  plasma,  the  (slightly  positive)  wall  potential  is: 

t  =^Mn[ri2  +  B]A(kTr/e)'] .  (12) 

t 

This  is  of  the  order  7^ ,  which  is  the  temperature  of 
the  secondary  electrons,  no  more  than  -  leV . 


=r,I2iT, -^}-T,8a[2kT„  -«&] .  (14) 

Thus,  for  Tr  <Tt  lP  and  Tc  > Tr  ^Jrespectively; 

%  =  Iexp^I2A-(T,-^rKC) 
h  e  kTt  ,•  (15) 


%-=i{2lfc(7,-7,J],  '  (16) 

J,  * 

where  jt  =  entc,l 4.  A  graph  of  this  is  shown  in  Fig. 
2.  Notice  the  sharp  increase  in  heat  loss  to  the  wall 
around  Ttgp .  In  fact,  the  heat  loss  grows  so  large  at 
TtJP  that  it  effectively  limits  the  electron  temperature 
7 i.tp  ■ 

Sheath  and  secondary  emission  effects  are  also 
important  when  considering  cross-field  electron 
transport  near  the  wall.  The  low-energy  secondary 
electrons  are  assumed  to  start  from  rest  at  the  wall  in 
crossed  electric  and  magnetic  fields.  By  calculating 
the  distance  traveled  downstream  by  their  guiding 
centers,  an  expression  for  the  total  near-wall  electron 
conductivity  is  determined: 


An  equation  for  electron  energy  lost  to  the  wall  and 
sheath  is  obtained  by  integrating  the  primary  and 
secondary  electron  energy  fluxes  across  Maxwellian 
distributions: 

=  f  »j/,(w)(f  m,w2)t/?w 
r  ,  ,  .  (13) 

-J  w\ftlwy±mtw*)d-  w’ 

where  n\  is  the  primary  electron  velocity  component 
normal  to  the  wall  and  where  w\  is  the  velocity  of 
the  secondary  electrons  assumed  to  be  Maxwellian  at 
temperature  Tv  .  The  integral  yields:  • 


Above,  8^  is  the  effective  secondary  emission  yield, 
and  &  is  the  angle  of  incidence  of  the  magnetic  field 
line  with  the  wall. 

Boundary  Conditions 

The  boundary  conditions  for  the  quasi- 1 -D  electron 
equations  are  handled  by  directly  fixing  7,  at  the 
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cathode,  and  by  imposing  a  zero-slope  condition  on 
Tt  at  the  anode. 

HPHal!  is  capable  of  modeling  background  (chamber) 
pressure  at  the  downstream  boundary4  however,  the 
background  pressure  is  set  to  zero  for  this  study.  The 
boundary  conditions  at  the  injector  also  include  the 
introduction  of  particles  from  the  propellant  feed. 
These  particles  are  placed  randomly  within  the 
injector  region,  and  take  random  trajectories 
corresponding  to  a  half-range  Maxwellian  at  a 
temperature  of  1000  degK. 


Numerical  Method 

The  governing  equations  are  solved  time-accurately  by 
separating  the  slow  time  scale  (ion  and  neutral)  motion 
from  the  fast  time  scale  (electron)  motion,  and 
iterating  successively.  Individual  ion  and  neutral 
aioms  are  simulated  using  a  Particle-ln-Cell  method. 
The  electron  motion  is  modeled  as  a  fluid  continuum 
with  the  differential  equations  derived  above  and 
solved  using  the  usual  methods. 

Runs  presented  in  this  paper  use  a  47-by-22  structured 
nonuniform  grid  encompassing  the  acceleration 
channel  and  approximately  4cm  of  the  plume. 
Rotational  symmetry  is  assumed.  Therefore,  only  a 
meridional  section  of  the  Hall  thruster  acceleration 
zone  is  modeled.  Grid  spacing  is  determined  by  the 
timestep  of  the  simulation.  U  is  set  not  to  exceed  the 
maximum  distance  traveled  by  an  ion  particle  in  one 
timestep. 

The  magnetic  field  is  generated  as  a  pre-process  by 
specifying  the  thruster  geometry,  assuming  infinite 
permeability  of  the  iron  poles,  and  solving  Laplace’s 
equation  on  the  regions  exterior  to  the  poles,  the 
coils  are  assumed  to  be  perfect  solenoids,  so  the 
problem  reduces  to  that  of  potential  flow,  with  each 
pole  piece  set  to  a  given  magnetic  potential. 

The  motion  of  heavy  particles  is  slow  compared  to  the 
electrons.  For  computational  efficiency,  die  electrons 
and  heavy  particles  are  moved  on  different  timesteps. 
Fig.  3  shows  the  sequence.  Since  Tt  does  not  vary 
along  magnetic  field  lines  (T,  =TtU)).it  is  possible 
to  reduce  the  electron  energy  equation  to  a  quasi-one- 
dimensiona!  form.  After  some  manipulation,  a  one¬ 
dimensional  nonlinear  differential  equation  is  derived 
for  Tr  as  a  function  of  X.  Hie  solution  of  this 
equation  is  accomplished  by  using  a  modified  Forward 
Time  Centered  Space  (FTCS)  method.1'  The  timestep 


J 

_ 3 

f 

5- 10'" 3 

Move  Ions  and  Neutrals  (PIC) 

v 

Ionize 

v 

Handle  Boundary  Conditions 

A/  =  5-10  s 


Output  Results 


Fig.  3.  Execution  sequence  of  the  numerical 
simulation. 


is  5-10'"  seconds,  based  on  successive  reduction  to 
the  case  where  the  solution  is  stable  and  unchanging 
for  smaller  timesteps.  Once  Te  and  d*  are  known  on 
the  domain.  6  is  found  using  Eqn.  (1). 


The  electric  field  and  electron  temperature, 
determined  by  integrating  the  electron  equations,  is 
then  used  in  a  PIC  method  for  ions  and  neutrals.  Ion 
positions  and  velocities  are  updated.  Also,  the 
densities  are  adjusted  appropriately  based  on 
computed  local  bulk  ionization  rate.  This  sequence 
repeats  as  shown  in  Fig.  3. 

Since  the  method  is  time-accurate,  the  simulation  will 
not,  in  general,  converge  to  a  steady  state  solution 
because  of  plasma  fluctuations.  Nevertheless,  a 
solution  is  considered  complete  when  the  fluctuations 
reach  a  regular  frequency  and  amplitude,  and  have 
repeated  many  periods.  For  this  paper,  results  are 
averaged  over  0.5ms .  One  case  takes  approximately 
5  hours  to  converge  on  a  Pentium  II  Xeon-class 
personal  computer. 


Results  and  Discussion 

Several  cases  were  run  to  examine  the  effect  of 
varying  the  secondary  electron  emission  of  the 
insulator  wall.  These  cases  are  shown  in  Table  1  in 
order  of  decreasing  secondary  electron  emission 
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Table  1.  Performance  results  from  the  numerical  simulation  for  various  secondary  electron  emission  coefficients 

of  the  wall  material. 


limns 

A 

N/A 

0.144 

0.126 

0.114 

0.104  ' 

0.096 

Tr.„<eV) 

N/A 

16.0 

20.0 

24.0 

28.0 

32.0 

itV) 

N/A 

15.5  ' 

19.03 

22.2 

25.0 

27.4 

tM ) 

123 

2.27 

2.40 

2.47 

2.52 

153 

4(A) 

1.61 

1.45 

1.53 

1-58 

1.60 

1.61 

FimN) 

37.8 

37.2 

39.3 

40.5 

41.1 

41.4 

n. 

0.938 

0.844 

0.891 

0.920 

0.932 

0.938 

0.722 

0.639 

0.638 

0.640 

0.635 

0.636 

1, 

0.673 

0.805 

0.806 

0.803 

0.806 

MEM 

1 

0.456 

0.434 

0.458 

0.473 

0.477 

0.483 

F‘* 

N/A 

0.052 

0.064 

0.074 

0.080 

0.086 

coefficient.  A.  from  (5).  The  exponential  component. 
B .  of  the  secondary  electron  emission  model  was 
constant  at  0.576. 

Geometry  is  that  of  an  SPT-70.  Operational 
parameters  are  Vj=30QV  and  m  =  2.34mgls .  The 
first  case  in  Table  1  is  experimental,  and  was  run  at 
AFRL.16  Some  experimental  values  are  marked 
"N/A"  in  Table  1  because  they  were  not  accurately 
measured.  The  remaining  independent  experimental 
measurands.  /„ .  fk.  and  F  have  estimated  accuracy 
of  2 9r.  10%.  and  5%,  respectively.6 

The  profile  of  electron  temperature  for  Case  1  is 
shown  in  Fig.  4.  The  electric  field  is  strongest  near 
the  exit  of  the  channel.  This  causes  Ohmic  heating  of 
the  electrons  in  this  region.  Their  energy  decreases 
closer  to  the  anode  because  of  inelastic  collision 
losses,  and  because  of  wall  losses.  The  contours  of 
electron  temperature  for  Cases  2-5  are  nearly  identical 
to  Case  1.  except  the  magnitude  increases  with 
decreasing  secondary  electron  emission  coefficient. 
The  peak  temperature.  Trjm ,  always  occurs  at  the 
exit,  and  is  given  for  each  case  in  Table  1 . 

The  increase  in  T„  with  decreasing  secondary  electron 
emission  coefficient  can  be  understood  in  the 
following  wa^*)  Assuming  wail  interactions  are  the 
primary  mechanism  of  electron  energy  loss.  T,  can  be 
expected  to  grow  until  limited  by  wail  losses.  As  can 
be  seen  in  Fig.  2.  the  heat  flux  to  the  walls  increases 
rapidly  as  Tr  approaches  Tr  tf .  Thus.  Te  „  becomes 
the  limit  of  electron  temperature.  Now,  from  (8)  and 
(II).  decreasing  the  secondary  emission  coefficient  of 


the  wall  material  increases  and.  therefore. 

* 

Efficiencies  shown  in  Table  1  are  for  the  -discharge 
only,  excluding  cathode  flow.  Breakdown  of  thrust 
efficiency  may  be  given  by 17  = BJU1,  •  where: 

em 

fc-f  (18) 

*a 

Above,  4  is  the  ion  beam  current.  /„  is  the  discharge 


Level  Tt(eV) 


Fig.  4.  Electron  temperature  contours  in  the  SPT-70 
discharge  channel  computed  by  HPHalL  Case  1. 
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Fig.  5.  Thrust  efficiency  versus  secondary  electron 
emission  coefficient.  For  Cases  1-5. 

current.  Vh  is  the  mean  ion  beam  energy  in  Volts,  and 
Vj  is  the  discharge  Voltage. 

As  secondary  electron  emission  coefficient  is 
decreased,  electron  temperature  increases,  as  does  die 
ionization  efficiency  of  the  plasma.  This  is  reflected 
in  the  growth  of  utilization  efficiency.  Tjr ,  with  T(  m 
for  Cases  1-5.  Fig.  5  shows  a  plot  of  total  efficiency 
versus  secondary  electron  emission  coefficient  A,  for 
Cases  1-5.  This  is  an  important  result,  because  it 
illuminates  a  predictable  connection  between 
secondary  electron  emission  and  thruster  performance. 

Another  noticeable  result  shown  in  Table  1  is  the 
increase  in  fraction  of  double  ions  in  the  beam, 
F"  -  IhMit  <  with  decreasing  secondary 

electron  emission  coefficient  A .  This  can  be 
attributed  to  the  cross  section  of  double  ionization  by 
electron  impact,  which  has  a  threshold  just  above  20 
eV.  The  small  changes  in  F"  observed  here  would 
have  a  negligible  effect  on  thruster  efficiency,  but  may 
be  an  important  consideration  for  studying  the 
interaction  of  Hall  thruster  beams  with  real  spacecraft 
components. 

Summary 

As  the  secondary  electron  emission  coefficient  was 
decreased  by  Is#:  HPHall  predicted  a  steady  increase 
in  Hall  thruster  efficiency  from  0.434  to  0.483.  This 
indicates  a  relationship  between  secondary  electron 
emission  of  Hall  thnister  insulator  wall  material  and 
thruster  performance.  It  suggests  that  secondary 
electron  emission  be  one  of  the  parameters  considered 
during  the  Hall  thruster  design  phase,  and  that 
materials  with  low  secondary  electron  emission  should 
be  favored. 
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